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Abstract: We argued in arXiv:1408.0624 that the quartic scalar held in AdS has 
features that could be instructive for answering the gravitational stability question of 
AdS. Indeed, the conserved charges identified there have recently been observed in the 
full gravity theory as well. In this paper, we continue our investigation of the scalar 
held in AdS and provide evidence that in the Two-Time Formalism (TTF), even for 
initial conditions that are far from quasi-periodicity, the energy in the higher modes 
at late times is exponentially suppressed in the mode number. Based on this and 
some related observations, we argue that there is no thermalization in the scalar TTF 
model within time-scales that go as ~ 1/e 2 , where e measures the initial amplitude 
(with only low-lying modes excited). It is tempting to speculate that the result holds 
also for AdS collapse. 
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1 Overview 

The question of whether AdS space is stable [1, 2] against turbulent thermalization 
and the formation of black holes under generic (non-linear) perturbations has received 
much attention recently. AdS space with conventional boundary conditions is like a 
box, and therefore perturbations that were weak to begin with can reflect multiple 
times from the boundary, potentially resulting in sufficient localization of energy to 
create black holes. Aside from the fact that black hole formation is a question of 
fundamental interest in (quantum) gravity, this problem acquires another interesting 
facet via the AdS/CFT correspondence: it captures the physics of thermalization in 
strongly coupled quantum field theories. 

At the moment however, it is fair to say that the evidence for and against the 
instability of AdS when excited by low-lying, low-amplitude modes is mixed [3-11]. 
In an effort to (partially) clarify this situation, in this paper we will make some 
comments about two loosely inter-related questions: 

• Does “most” initial data lead to thermalization? 

• Can one argue that within a time-scale of order 0(l/e 2 ), where e captures 
the amplitude of the initial pertubration, thermalization does (not?) happen? 
This is an interesting question because the statement of [1] is that black hole 
formation happens in AdS within this time scale. 

We will ask these questions, which are inspired by gravitational (in)stability in AdS, 
in the context of a simpler problem: a self-interacting 0 4 scalar field in AdS. The 
works of [6-10] suggest that these systems have close similarities, so we believe that 
this effort will be instructive and worthwhile. 

One of our main tools will be the Two-Time Formalism (TTF) developed in [4] 
(we will describe this approach in section 2). We will argue why this approach has 
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various advantages, and why we believe it captures the essential physics of resonances 
in the full (ie., non-TTF) model. But we emphasize that this will shed light on the 
instability question only if the instability, if it exits, is caused by resonances (which 
seems plausible to us). If the instability is caused by some other (possibly longer 
time-scale) dynamics, TTF in the leading order can miss that physics. But we expect 
that physics in the 0( 1/e 2 * ) time-scale should be captured by TTF. 

Furthermore, for concreteness, we will take the following as the definition of 
the absence of thermalization : the presence of exponentially distributed energies in 
the higher modes, as a function of the mode number 1 . That is, if the system has 
Aj ~ e ~at late times for some positive (3 we will say that it does not thermalize 
(at least for a very long time). Loosely, one could also adopt a definition that the 
system has instability towards thermalization if the late-time behavior of the j’th 
mode goes as Aj ~ j~ a , where a is a positive quantity — it is possible however that 
this is not a necessary nor a sufficient condition [8], and we will not use this in our 
paper. 

Within the context of these three limitations (namely, working with the scalar 
field and in the TTF approach and within a particular definition of thermalization) 
our results imply the following answers for the two questions (combined into one): 

• Initial data with only low lying modes do not lead to thermalization for the 
quartic scalar field in the TTF formalism within a time-scale of 0( 1/e 2 ). This 
suggests that if at all there is thermalization in the full theory, it should be 
coming from non-resonant transfer of energy. 

Together, we believe that these observations present fairly strong evidence that 
thermalization (as defined above) does not happen for initial value data which have 
only the low-lying modes excited. Our results, as already emphasized, are for the 
0 4 -scalar: but we believe similar statements apply for AdS gravity as well. We make 
various further comments of varying degrees of technicaility in later sections. 

This result might seem at odds with the numerical results for the Gaussian initial 
value profile discussed in [1], But it has been pointed out in [8] that the spectrum 
of the Gaussian profile in AdS, in fact has energy ditributed in the higher modes 
as a power law to begin with. So it is not a suitable initial initial profile to resolve 
between the presence or absence of collapse using our criterion for the absence of 
thermalization, which requires exponential suppression in the higher modes. Clearly 

1 Note that the definition of thermalization is somewhat ambiguous. We are adopting this as 

a sufficient but not necessary condition for the absence of thermalization as we will make more 

precise at the beginning of Section 3. One source of ambiguity is that our system is classical and 

suffers from a UV catastrophe: so once the system has fully thermalized, the average energy per 
state would be zero, if we don’t truncate it. In particular, the distribution of energies should not 
be compared to a canonical ensemble distribution, rather it should be thought of as capturing the 
efficiency of energy transfer to higher modes. 
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more work is required to clarify whether this is an acceptable initial profile or not 
for resolving this particular question. 

For completeness, lets also state that our results are still not quite conclusive. 
Apart from the points emphasized above, there is also the perverse possibility that 
collapse happens, but not due to resonances - but note however that the time-scale 


for this will be bigger than ~ 1/e 2 . 

2 TTF Formalism 

The action for the scalar field theory is given by 

S = J d * xV=g + V (</,)) (2.1) 

where the potential is given by 

VW = j/ ( 2 . 2 ) 

The metric for the space is given by 

ds 2 = sec 2 x (— dt 2 + dx 2 + sin 2 x dQ 2 ) (2.3) 

The equations of motion for the scalar field are given by 

0(2,0) + Dg 0 = 0(2,0) _ 0(0,2)-2- m = -A 0 3 

sin x cos x 6 cos 2 x 

where the D s represents the spatial Laplacian operator. This operator has an eigen¬ 
function basis given by 

n s ej(x) = lu 2 ej(x) (2.5) 

e/x) = 4^ ^ + + ^ cos 3 x 2 F\ + 3; sin 2 x^ (2.6) 

u] = (2j + 3) 2 j = 0,1,2,... (2.7) 

The inner product in this basis is defined as 

( f,g)= [ dx tan 2 x f(x) g(x) (2.8) 


In the Two-Time Framework (TTF), we have the slow-moving time defined as r = 
e 2 t, which requires the time derivatives to be redihned as dt —> dt + e 2 d T . The scalar 
field is written as an expansion in the small-parameter e as 

(j) = e 0(i) (t, r, x) + e 3 0 (3 ) (t, r, x) + £>(e 5 ) (2.9) 
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Note that the ratio between the slow and fast times (r and t) also controls the overall 
scale of the amplitude. Putting this expansion in the scalar field equations of motion 
eq.(2.4) we get 

2 

order e: <9 2 </>m(f, t, x) - d%<f>m(t, r, x) - : - d x <f>m(t, r, x) = 0 (2.10) 

Sill x cos x 

order e 3 : <9 2 0 (3) (t, r, x) + 2d t d T (() (1) (t, r, x) - <9 2 0 (3) (t, r, x) 

2 „ . A 


(<9x0(3 )(t,T,x) = 


Sill X COS X 

The order e equation has the general real solution 


6 cos 2 x 


0(1 )(t,T,x) (2.11) 


0(1) {t, T,x) = ^2 ( A j( T ) e iU}jt + A j(,T)e lUljt ) 6j{x) 

3=0 


( 2 . 12 ) 


Note that the introduction of the slow times gives an extra variable that we can tune 
- we will use this at order e 3 to cancel of the resonant terms. The equations that 
accomplish this are called the TTF equations. Substituting the above first order 
results into the order e 3 equations we get 

OO 

d?(j) {3) (t,T,x) - 2i (a r A i (T)e“* Wj ' t - d T Aj(T)e Wit ) e^x) + □ s 0( 3 )(i,T,x) 

A 


6 cos 2 x 


k =0 

oo 

E 

j,k,l =0 . 


+ Aj(T)e iu 2 (A k (r)e~ iuJkt + A k (r)e iUkt ) 


(Ai(r)e lUlt + Ai(r)e lUlt ) ej(x)e k (x)ei(x) 

Projecting on the basis solutions give 

ej(x), [<9 2 + (u|]0 (3) (t,r,x)) - 2 iuj [d T Aj(T)e~ luljt - d T Aj(T)e luJjt ] 


(2.13) 


A 


E c 


jklm 


k,l,m =0 


\A k {r)e~ iUkt + A k (r)e^ kt ] [M^e^ + A z (r)e^] 


[Am^e-^ + A m {r )e iuJmt ] 


where 


C'jklm 


t/2 


dx tan 2 (x) sec 2 (x) ej(x)e k (x)ei(x)e m (x) 


(2.14) 


(2.15) 


By direct computation (using properties of Jacobi polynomials - which are an alter¬ 
nate way to describe the basis functions, see Appendix A), one can show that the 
necessary and sufficient condition for resonances is 


Wj + U3 m — u k + UJi. 


(2.16) 
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The absence of other combinations for the resonances for the scalar theory was rec¬ 
ognized and used in [ 6 ] (see footnote 3 of [7] for a simple proof). They are also 
absent in the gravity case, but the computation required to show this in that case is 
substantially more lengthy [5]. The close parallel between the structure of the reso¬ 
nances in the two cases is evidently one of the reasons why they exhibit similarities 
in their thermalization dynamics [ 6 ]. 

In any event, at this stage we have the freedom to choose the Aj(t) as mentioned 
above so that the resonances on both sides are cancelled. This is accomplished by 
solving the Aj according to 

A °° 

— 2iu>j d r Aj = — — C jk i m A k AiA m (2-17) 

k,l,m =0 

and its complex conjugate. These are the TTF equations that we will use extensively 
in the next section. Once the resonances are cancelled, the coupling to the higher 
modes is expected to be weak and we believe it is unlikely that there will be efficient 
channels for thermalization: but this is a prejudice, and possibly far from proof. In 
any event, we can systematically solve for 0 ( 3 ) (t, r, x) at this stage if we wish, without 
being bothered by resonances. 

Note that the simplicity of the quartic scalar arises from the fact that the C i]k i 
have a (relatively) simple expression. We will comment more about this in Appendix 
A. 

Before concluding this section we quote some pertinent results from [ 6 ] for our 
scalar TTF system. Firstly, we can get the TTF equations using an effective La- 
grangian 

L TT f = - BiBi) -\-^C ijk iB i (T)B j (T)B k (T)Bi( y T), (2.18) 

i 

where summation in the interaction term is over i,j,k,l such that cU; + ujj — u> k — 0 Ji = 
0. In writing the expression in this form, we have done an appropriate scaling of 
each mode by uj k and A for easy comparison with the notation of [ 6 ]: B k are the 
rescaled modes. The system has a dilatation symmetry: B k (r) —>■ eB k (^r). So if 
thermalization happens in the TTF theory it should scale inversely as the square of 
the amplitude: the assumption that TTF theory captures the relevant physics is the 
assumption that the system has such a scaling regime. 

However, the system has the following conserved charges [ 6 ] arising from a cor¬ 
responding set of symmetries: 

Qo = ^B k B k , symmetry : B k ->• e* e B kl (2.19) 

Qi = ^ kB k B k , symmetry : B k ->■ e lke B k , (2.20) 

E = ^2 Cpfc«5 i (r)i? i (r)H fc (r)H z (r), Symmetry : t ->■ t + a. (2.21) 

uJi -\-Uj — cjfc —to i =0 
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Various pieces of evidence indicating that the evolution of the quartic scalar in AdS 
has some close connections to collapse in AdS gravity were presented in [6]. The 
above conserved charges were identified for the full gravity system in [7] (see also 
[ 8 ])- 

3 Results 

In this section, we will study various aspects of the TTF equations for the quartic 
scalar in some detail. As mentioned in the introduction, we will take the exponential 
decay of Aj(r) with j as an indication that thermalization is suppressed. In [8] 
some arguments were made that a power law Aj ~ j~ a for positive a is indicative 
of thermalization/black hole formation. We will make this somewhat more precise 
as follows. The basic object that is taken as an indicator of collapse in [1, 4] is 
the quantity |II(t,0)| 2 , the unbounded growth of whose profile is taken as the onset 
of collapse. The analogue of this quantity in our scalar TTF case can be taken as 
|0(i)(f, 0)| J (compare Figure 1(A) and the accompanying discussion in [6] to Figure 
3 in [4]). At this point, using (2.12), (2.6) and (2.7) we can see that this quantity 
can be estimated and bounded via 

|0d)(t,o)T ~ 15T 2/ tf ~ (3-D 

Now, it is evident that when Aj ~ e~^ a the last quantity is finite and therefore 
the LHS can never diverge, which is what we set out to show. This indicates that 
expoential suppression of higher modes is a sufficient condition for absence of ther¬ 
malization. Note however that we are silent about what constitutes thermalization 
at the level of modes - fortunately, we will never need a precise definition of that for 
the purposes of this paper. 

The TTF theory has quasi-periodic solutions (see [4] for a discussion of analogous 
solutions in the gravity system) of the form 

Aj(r) = aj exp(—z/3j-r), where f3j = /3 0 + j(/3 1 - /3 0 ). (3.2) 

One can choose ao,a± (or f3o,/3±) and determine the rest of the aj via the TTF 
equations (2.17) 2 , if one truncates the system at some j = j max and demand stability 
of the solution against variation in j max ■ We have done this, and the resulting modes 
decay with the mode number j as ~ cxp (r c T f or some positive c, see Figure (1). This 
is obviously consistent with our definition of (non-)thermalization. In [4] the jmax 
was taken to be ~ 50, in our case we are able to go up to j max = 150. 

2 For some initial conditions we see more than one quasiperiodic solution. 
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Figure 1: The log-plot of jotj vs. j for the quasi-periodic solutions. The linear fit 
is indicative of exponential suppression of Aj with j. 


Logll^ll 



(a) Representative plots of “small” (b) Representative plots of “large” 
perturbations around quasi-periodic perturbations around quasi-periodic 
solutions. solutions. 

Figure 2: Evolution plots of perturbations around quasi-periodic solutions. 



If we perturb a quasi periodic solution we expect to get oscillations of the A/s 
around ay. See Figure (2a) for solutions where the initial value of the Aj are close 3 to 
their quasiperiodic values. If on the other hand, the initial Aj values are sufficiently 
far from their quasiperiodic values, we expect that the solutions transition to chaos. 
This expectation is qualitatively verified in Figure (2b) where we launch the Aj 
far away from quasi-periodicity. In what follows we will show that even in these far- 
frorn quasiperiodic solutions, the maximum value attained by the Aj as we evolve the 
solution is exponentially suppressed in j. This is an indication that energy transfer 
to the higher modes is suppressed even in these solutions - if this behavior holds also 

3 In order to make these statements precise, we will need a notion of closeness between solutions 
in terms of modes. A convenient way to define a dimensionless measure of the “distance” between 
two solutions (say 1 and 2) is to consider 


Aio — 


E.M^. 


( 1 ) A ( 2 )* 


VE 




Vej4 2) i 2 


A 12 ~ 1 is close. The summation is only up to mode number j max . 


(3.3) 
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(a) i,j, k, l = 1, 2, n, n + 1. (b) i, j, k,l = n, 2 n, 3n, 4n. 

Figure 3: Plot of Cijki together with an inverse linear fit. 


in gravity, it conld be an indication that these solutions generically do not collapse. 

One of the ways in which one might try to understand the efficiency of energy 
transfer to higher modes is by studying the coefficients 4 C^ki which signify the cou¬ 
pling between the modes. To understand the behaviour of TTF equations at large 
j, we look at various kinds of limits we may consider for as the i,j,k,l are 
sent to oo. One is a simple scaling of indices, i,j,k,l —> ai,aj,ak,al. By fitting 
the plot (see Figure 3b), we see that in this case C t] ki goes as O(-) as a —» oo. 
Another case is where we keep two modes fixed and take another two to infinity: 
i ~ j ~ approximately fixed, but with k ~ l ~ a and we take a —> oo. We fold that 
they also have a O(-) fall off. It is important to note that because of the resonance 
condition, these are the only possible couplings available for a high mode - one can¬ 
not (for example) hold three indices small while sending the forth one to infinity. So 
progressively higher modes are weakly coupled, both to each other as well as to the 
low-lying modes. 

Finally, we consider the evolution of the modes when we launch the system both 
near and far from quasi-periodic initial conditions. The way we do this is by calcu¬ 
lating the coefficients C^ki analytically (see Appendix A for some comments on this) 
and then integrating the resulting TTF equations numerically for the various initial 
conditions. In all cases we plot maximum value of A 3 that is attained during the 
entire period of evolution against j, and we fold that this MaxfA^r)] exponentially 
decays with j for all initial data. We see an exponential decay with respect to j, 
not just for solutions close to quasi-periodic solutions, but also for those that are 
far from it: see Figure (4). This is true even though for some initial conditions we 
see a approximately power law decay of modes up to some intermediate frequency. 
These statements can be verified using the norm (3.3) with the understanding that 
the summation over j has to be restricted to be above some appropriately chosen 
jmin (and of course below j max ) when we are talking about high modes. 

4 Note that the coefficients C^ki can be determined via (2.15) analytically, but using Mathemat- 
ica. Some comments on this are given in Appendix A. 
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Figure 4: Plot of log[Max[A n (r)]] for j G [0,150] and a linear fit. The fit has been 
done in the region j G [40,150]. The last figure corresponds to quasi-periodic initial 
data. 
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A Comments on Cijki and Jacobi Ploynomials 


The determination of the C^ki is in principle straightforward by direct evaluation of 
(2.15). This is an analytically tractable problem because the basis functions ej(x) 
can be written in terms of Jacobi polnomials as 


i( x ) = 4 


(j + i)U + 2 ) r(j + i)r(3/2) 
7T r(j + 3/2) 


cos 3 x P 


( 1 / 2 , 3 / 2 ) 


(cos2x). (A.l) 
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Jacobi polynomials are (orthogonal) polynomials in their arguments and therefore in 
onr case they merely involve only (a finite number of) powers of sinusoids 5 . There¬ 
fore the integral for Cijki, which is in the range [0, 7 t/ 2] can, again in principle, be 
straightforwardly evaluated. It turns out that the result can be expressed in terms 
of finite sums of finite products of Gamma functions and such, but simplifying them 
on Mathematica becomes time-consuming. One could in principle try to simplify the 
expressions manually, but we have adopted a more pragmatic approach: we evaluate 
the integrals analytically on Mathematica by re-expressing the powers of sinusoids 
in terms of product formulas. Since the integrals are over [0, 7 t/ 2] this makes them 
substantially less intensive as far as time requirements are considered. This way we 
are able to algorithmize the (analytic) computation of C^ki on Mathematica, after 
which we use them in the TTF equations to do our numerical evolutions. 
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